Cortical lamina-dependent blood volume changes in human brain at 7T

Cortical layer-dependent high (sub-millimeter) resolution functional magnetic resonance imaging (fMRI) in human or animal brain can be used to address questions regarding the functioning of cortical circuits, such as the effect of different afferent and efferent connectivities on activity in specific cortical layers. The sensitivity of gradient echo (GE) blood oxygenation level-dependent (BOLD) responses to large draining veins reduces its local specificity and can render the interpretation of the underlying laminar neural activity impossible. The application of the more spatially specific cerebral blood volume (CBV)-based fMRI in humans has been hindered by the low sensitivity of the noninvasive modalities available. Here, a vascular space occupancy (VASO) variant, adapted for use at high field, is further optimized to capture layer-dependent activity changes in human motor cortex at sub-millimeter resolution. Acquired activation maps and cortical profiles show that the VASO signal peaks in gray matter at 0.8–1.6 mm depth, and deeper compared to the superficial and vein-dominated GE-BOLD responses. Validation of the VASO signal change versus well-established iron-oxide contrast agent based fMRI methods in animals showed the same cortical profiles of CBV change, after normalization for lamina-dependent baseline CBV. In order to evaluate its potential of revealing small lamina-dependent signal differences due to modulations of the input-output characteristics, layer-dependent VASO responses were investigated in the ipsilateral hemisphere during unilateral finger tapping. Positive activation in ipsilateral primary motor cortex and negative activation in ipsilateral primary sensory cortex were observed. This feature is only visible in high-resolution fMRI where opposing sides of a sulcus can be investigated independently because of a lack of partial volume effects. Based on the results presented here, we conclude that VASO offers good reproducibility, high sensitivity and lower sensitivity than GE-BOLD to changes in larger vessels, making it a valuable tool for layer-dependent fMRI studies in humans

Vascular autorescaling of fMRI (VasA fMRI) improves sensitivity of population studies : A pilot study

The blood oxygenation level-dependent (BOLD) signal is widely used for functional magnetic resonance imaging (fMRI) of brain function in health and disease. The statistical power of fMRI group studies is significantly hampered by high inter-subject variance due to differences in baseline vascular physiology. Several methods have been proposed to account for physiological vascularization differences between subjects and hence improve the sensitivity in group studies. However, these methods require the acquisition of additional reference scans (such as a full resting-state fMRI session or ASL-based calibrated BOLD). We present a vascular autorescaling (VasA) method, which does not require any additional reference scans. VasA is based on the observation that slow oscillations (<0.1Hz) in arterial blood CO2 levels occur naturally due to changes in respiration patterns. These oscillations yield fMRI signal changes whose amplitudes reflect the blood oxygenation levels and underlying local vascularization and vascular responsivity. VasA estimates proxies of the amplitude of these CO2-driven oscillations directly from the residuals of task-related fMRI data without the need for reference scans. The estimates are used to scale the amplitude of task-related fMRI responses, to account for vascular differences. The VasA maps compared well to cerebrovascular reactivity (CVR) maps and cerebral blood volume maps based on vascular space occupancy (VASO) measurements in four volunteers, speaking to the physiological vascular basis of VasA. VasA was validated in a wide variety of tasks in 138 volunteers. VasA increased t-scores by up to 30% in specific brain areas such as the visual cortex. The number of activated voxels was increased by up to 200% in brain areas such as the orbital frontal cortex while still controlling the nominal false-positive rate. VasA fMRI outperformed previously proposed rescaling approaches based on resting-state fMRI data and can be readily applied to any task-related fMRI data set, even retrospectively

Layer-specific activation of sensory input and predictive feedback in the human primary somatosensory cortex

When humans perceive a sensation, their brains integrate inputs from sensory receptors and process them based on their expectations. The mechanisms of this predictive coding in the human somatosensory system are not fully understood. We fill a basic gap in our understanding of the predictive processing of somatosensation by examining the layer-specific activity in sensory input and predictive feedback in the human primary somatosensory cortex (S1). We acquired submillimeter functional magnetic resonance imaging data at 7T (n = 10) during a task of perceived, predictable, and unpredictable touching sequences. We demonstrate that the sensory input from thalamic projects preferentially activates the middle layer, while the superficial and deep layers in S1 are more engaged for cortico-cortical predictive feedback input. These findings are pivotal to understanding the mechanisms of tactile prediction processing in the human somatosensory cortex

Investigation of the neurovascular coupling in positive and negative BOLD responses in human brain at 7T

Decreases in stimulus-dependent blood oxygenation level dependent (BOLD) signal and their underlying neurovascular origins have recently gained considerable interest. In this study a multi-echo, BOLD-corrected vascular space occupancy (VASO) functional magnetic resonance imaging (fMRI) technique was used to investigate neurovascular responses during stimuli that elicit positive and negative BOLD responses in human brain at 7 T. Stimulus-induced BOLD, cerebral blood volume (CBV), and cerebral blood flow (CBF) changes were measured and analyzed in ‘arterial’ and ‘venous’ blood compartments in macro- and microvasculature. We found that the overall interplay of mean CBV, CBF and BOLD responses is similar for tasks inducing positive and negative BOLD responses. Some aspects of the neurovascular coupling however, such as the temporal response, cortical depth dependence, and the weighting between ‘arterial’ and ‘venous’ contributions, are significantly different for the different task conditions. Namely, while for excitatory tasks the BOLD response peaks at the cortical surface, and the CBV change is similar in cortex and pial vasculature, inhibitory tasks are associated with a maximum negative BOLD response in deeper layers, with CBV showing strong constriction of surface arteries and a faster return to baseline. The different interplays of CBV, CBF and BOLD during excitatory and inhibitory responses suggests different underlying hemodynamic mechanisms

Functional cerebral blood volume mapping with simultaneous multi-slice acquisition

The aim of this study is to overcome the current limits of brain coverage available with multi-slice echo planar imaging (EPI) for vascular space occupancy (VASO) mapping. By incorporating simultaneous multi-slice (SMS) EPI image acquisition into slice-saturation slab-inversion VASO (SS-SI VASO), many more slices can be acquired for non-invasive functional measurements of blood volume responses. Blood-volume-weighted VASO and gradient echo blood oxygenation level-dependent (GE-BOLD) data were acquired in humans at 7 T with a 32-channel head coil. SMS-VASO was applied in three scenarios: A) high-resolution acquisition of spatially distant brain areas in the visuo-motor network (V1/V5/M1/S1); B) high-resolution acquisition of an imaging slab covering the entire M1/S1 hand regions; and C) low-resolution acquisition with near whole-brain coverage. The results show that the SMS-VASO sequence provided images enabling robust detection of blood volume changes in up to 20 slices with signal readout durations shorter than 150 ms. High-resolution application of SMS-VASO revealed improved specificity of VASO to GM tissue without contamination from large draining veins compared to GE-BOLD in the visual cortex and in the sensory-motor cortex. It is concluded that VASO fMRI with SMS-EPI allows obtaining a reasonable three-dimensional coverage not achievable with standard VASO during the short time period when blood magnetization is approximately nulled. Due to the increased brain coverage and better spatial specificity to GM tissue of VASO compared to GE-BOLD signal, the proposed method may play an important role in high-resolution human fMRI at 7 T

Sub-millimeter fMRI reveals multiple topographical digit representations that form action maps in human motor cortex

The human brain coordinates a wide variety of motor activities. On a large scale, the cortical motor system is topographically organized such that neighboring body parts are represented by neighboring brain areas. This homunculus-like somatotopic organization along the central sulcus has been observed using neuroimaging for large body parts such as the face, hands and feet. However, on a finer scale, invasive electrical stimulation studies show deviations from this somatotopic organization that suggest an organizing principle based on motor actions rather than body part moved. It has not been clear how the action-map organization principle of the motor cortex in the mesoscopic (sub-millimeter) regime integrates into a body map organization principle on a macroscopic scale (cm). Here we developed and applied advanced mesoscopic (sub-millimeter) fMRI and analysis methodology to non-invasively investigate the functional organization topography across columnar and laminar structures in humans. Compared to previous methods, in this study, we could capture locally specific blood volume changes across entire brain regions along the cortical curvature. We find that individual fingers have multiple mirrored representations in the primary motor cortex depending on the movements they are involved in. We find that individual digits have cortical representations up to 3 ​mm apart from each other arranged in a column-like fashion. These representations are differentially engaged depending on whether the digits’ muscles are used for different motor actions such as flexion movements, like grasping a ball or retraction movements like releasing a ball. This research provides a starting point for non-invasive investigation of mesoscale topography across layers and columns of the human cortex and bridges the gap between invasive electrophysiological investigations and large coverage non-invasive neuroimaging

High-resolution CBV-fMRI allows mapping of laminar activity and connectivity of cortical input and output in human M1

Layer-dependent fMRI allows measurements of information flow in cortical circuits, as afferent and efferent connections terminate in different cortical layers. However, it is unknown to what level human fMRI is specific and sensitive enough to reveal directional functional activity across layers. To answer this question, we developed acquisition and analysis methods for blood-oxygen-level-dependent (BOLD) and cerebral-blood-volume (CBV)-based laminar fMRI and used these to discriminate four different tasks in the human motor cortex (M1). In agreement with anatomical data from animal studies, we found evidence for somatosensory and premotor input in superficial layers of M1 and for cortico-spinal motor output in deep layers. Laminar resting-state fMRI showed directional functional connectivity of M1 with somatosensory and premotor areas. Our findings demonstrate that CBV-fMRI can be used to investigate cortical activity in humans with unprecedented detail, allowing investigations of information flow between brain regions and outperforming conventional BOLD results that are often buried under vascular biases

LayNii: a software suite for layer-fMRI

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain ‘layerification’ and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data